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In this paper, we develop the analysis of a two-dimensional magnetohydrodynamical configuration 
for an axially symmetric and rotating plasma (embedded in a dipole like magnetic field), modeling 
the structure of a thin accretion disk around a compact astrophysical object. Our study investigates 
the global profile of the disk plasma, in order to fix the conditions for the existence of a crystalline 
morphology and ring sequence, as outlined by the local analysis pursued in [T] [2]. In the linear 
regime, when the electromagnetic back-reaction of the plasma is small enough, we show the existence 
of an oscillating radial behavior for the flux surface function which very closely resembles the one 
outlined in the local model, apart from a radial modulation of the amplitude. In the opposite 
limit, corresponding to a dominant back-reaction in the magnetic structure over the field of central 
object, we can recognize the existence of a ring-like decomposition of the disk, according to the 
same modulation of the magnetic flux surface, and a smoother radial decay of the disk density, with 
respect to the linear case. In this extreme non-linear regime, the global model seems to predict a 
configuration very close to that of the local analysis, but here the thermostatic pressure, crucial for 
the equilibrium setting, is also radially modulated. 

Among the conditions requested for the validity of such a global model, the confinement of the 
radial coordinate within a given value sensitive to the disk temperature and to the mass of the 
central objet, stands; however, this condition corresponds to dealing with a thin disk configuration. 
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I. INTRODUCTION 

The understanding of the accretion process [3] char- 
acterizing compact astrophysical objects is regarded as 
a central issue when describing the triggering of high 
energy phenomena, such as the formation of jets. In- 
deed, it is commonly believed that the intense electro- 
magnetic emission in different bands from astrophysical 
sources, such as Gamma Ray Bursts (see for example [I]) 
or Active Galactic Nuclei (see for example [5]), should be 
the final result of the proper instabilities of the accretion 
mechanism. Theoretical attempts have been made over 
the years to explain the formation of jets from the accre- 
tion profile of compact and massive astrophysical systems 
(see, among others, [B] and [?])■ 

The accretion mechanism in astrophysical systems is, 
however, far from being settled, because many aspects of 
the angular momentum transport within the gravitating 
plasma are still open to scientific debate. The stellar disk 
configurations, on which the present analysis is focused, 
are described by a well-established paradigm, which, for 
the case of a thin disk profile, has the character of a stan- 
dard model pJj- This scenario satisfactorily reproduces 
the phenomenology fixed by the observational data, but 
at the price of postulating a viscoresistive magnetohydro- 
dynamic (MHD) approach, which is not supported by the 
microscopic features of the disk plasma. The possibility 
to reconcile such an effective description of the non ideal 
nature of the plasma with the microscopic estimations of 
the viscosity and resistivity coefficients is individuated in 
the appearance of strong turbulence behaviors inside the 
disk. The theoretical framework for describing the on- 
set of the plasma instabilities is offered by the so-called 
magnetorotational instability, introduced in [U [9] (see 

also mam])- 

Even if the standard model for the thin accretion disks 
is rather successful and widely accepted, it cannot es- 
cape the following two relevant criticisms, (i) The z- 
dependence of the configuration is somewhat frozen out 
from the equilibrium by an average procedure along the 
vertical direction of both the radial and azimuthal equa- 
tions, as originally proposed in [12j[l3]. (ii) The experi- 
ence acquired in plasma morphology from the laboratory 
activities, does not confirm the direct relation between 
the existence of a turbulent regime and the applicability 
of the viscoresistive MHD scenario [14] . 

In |15j it is shown how the configuration of a thin disk, 
treated in a two-dimensional ideal MHD scheme, reveals 
very different features with respect to the effective one- 
dimensional standard model. In fact, a new equilibrium 
profile, corresponding to a strong magnetization of the 
disk, is determined by virtue of the confinement induced 
by the Lorentz force. The main significant property of 



this regime is the non-Keplerian character of the angular 
velocity of the disk. 

A more specific criticism to the standard model is pur- 
sued in p] [2] , where a fundamental description of the ro- 
tating disk plasma is addressed by postulating the ideal 
two-dimensional MHD and arguing that the accretion 
mechanism can be driven by intermittent ballooning in- 
stabilities which push the plasma toward the central ob- 
ject via its porosity near the X-points of the magnetic 
configuration |16| . Until now, the main achievement of 
this new point of view, is the demonstration that the 
local configuration around a fixed value of the radial co- 
ordinate, outlines the decomposition of the disk into a se- 
quence of rings, corresponding to adjacent and opposite 
current density filaments. This issue takes place only in 
the limit of a very strong electromagnetic back-reaction 
of the plasma (its existence in the limit of large values 
of the plasma j3 parameter has been investigated in |17|). 
but a periodic (crystalline) structure of the magnetic sur- 
face functions already emerges in the linear regime. 

In the present analysis, we fix the conditions necessary 
to recover the crystalline and ring like profile within a 
global two-dimensional MHD model. In other words, we 
quantify the role played by the radial envelope in deter- 
mining the amplitude scaling of the oscillating behavior 
that characterizes the magnetic configuration. 

In the linear case, we are able to sketch a clear pic- 
ture for the emergence of the crystalline profile from the 
disk configuration, establishing that the amplitude of the 
perturbed magnetic surface scales as ^/r, while the mass 
density on the equatorial plane decays as r~ 3 (r being the 
radial coordinate). In order to obtain this scheme from 
the generic axially symmetric MHD scenario, we need to 
make some significant assumptions, which, however, can 
be mainly summarized by the thin nature of the disk, the 
short scale of the perturbations in the plasma, and by the 
confinement of the configuration well-inside a given ra- 
dial region, determined by the model parameters (such as 
the mass of the central object and the disk temperature). 

When the electromagnetic back-reaction within the 
disk plasma induces a dominant magnetic field (with re- 
spect to the one provided by the central object), we re- 
cover the same ring like decomposition described in the 
local model [2]. The radial modulation of the amplitude 
of the oscillations of the magnetic surfaces is the same as 
in the linear case, but now the mass density of the disk 
decays slower along the radial direction, i.e., like r -1 / 2 
instead of r~ 3 . Furthermore, in the non linear case, the 
pressure assumes a significant role in fixing the equilib- 
rium, and its oscillation amplitude is also modulated as 
r -1 , while retaining the same structure as in the local 
model. 

The main achievement of the present analysis is to 
show that the local features outlined in [T] for the steady 
configuration of a thin disk have a global character too 
and therefore they may concern fundamental astrophys- 
ical processes, taking place within the disk plasma, such 
as the realization of material jets [T5J |T5] . Furthermore, 
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we provide a precise constraint about the internal radial 
region in which the crystal profile of the magnetic field 
can take place, identifying a parameter discriminating for 
specific astrophysical sources. 

The paper is organized as follows. In Sec. [TTJ we give 
the motivation for a reformulation of the accretion prob- 
lem toward a more ideal nature of the gravitating plasma. 
In Sec. |III| we review the fundamental equations of 
axially-symmetric two-dimensional MHD describing an 
accretion disk embedded in the gravitational and mag- 
netic fields of a central object; furthermore we fix the 
relation between the angular frequency of the plasma 
and the flux function. In Sec. |IV[ we develop the per- 
turbation scheme for the considered problem, expand- 
ing the configuration equations up to the first order in 
agreement with the assumption of dealing with small- 
scale perturbations, and fix the conditions needed to re- 
cover the crystalline structure in the global case. In Sec. 
\V\ we assess the linear case when the back-reaction of 
the plasma is much smaller than the dipole-like field of 
the central object, and, in Sec. |VI| we discuss the phe- 
nomenological implications of the constructed model in 
order to outline the main features of the global plasma 
configuration. Then, in Sec. |VII| we analyze the opposite 
case, i.e., when the back-reaction of the plasma cannot 
be considered as a small perturbation but it is the domi- 
nant contribution in fixing the global profile. Concluding 
remarks follow in Sec. IVIIII 

II. NECESSITY FOR A REFORMULATION OF 
THE ACCRETION PICTURE 

In the Standard Model for the stellar accretion [31 [T2] , 
the configuration of an axisymmetric thin disk is deter- 
mined by the fluidodynamical equilibrium which takes 
place in the gravitational field of the compact accreting 
object (having mass M*), but performing an average pro- 
cedure along the vertical direction. 

The radial equilibrium states that the angular fre- 
quency of the disk takes the Keplerian profile ui(r) = 
cjk = \J GAY* /r 3 . Significant deviations from such a be- 
havior are expected only in advective dominated regimes. 

The vertical equilibrium in the isothermal disk of tem- 
perature T, fixes the exponential decay of the mass den- 
sity p over the equatorial plane value po{r), i.e., p/po = 
D(z 2 ) = exp{— z 2 /H 2 }, where we introduced the typical 
length H 2 = 2v 2 /u!^ estimating the half-depth of the disk 
(v s is the sound velocity on the equatorial plane, namely 
v 2 = 2K-qT jm u with Kb the Boltzmann constant and rti\ 
the ion mass). The azimuthal equilibrium describes the 
angular momentum transport across the disk, by virtue 
of a turbulent viscosity coefficient D, such that 

M d (L - L d ) = 3irT)uj K r 2 , (1) 

where L is the angular momentum per unit mass, L d is 
a fixed value and A Yd = —2wrHv r is the mass accretion 



rate, associated with the radial equatorial velocity v r < 
and the surface mass density E = J_ H pdz. Finally, 

the continuity equation implies that M d = const > (a 
discussion of these equations can be found, for example, 
in [3]). 

A. The viscoresistive puzzle 

The microscopic plasma structure accounts for a too 
small viscosity coefficient D arising in the disk to ex- 
plain the accretion rates observed in some astrophysical 
systems, such as X-ray binaries. In fact, the observed ac- 
cretion rates, evaluated by the increasing disk luminosity 
id ~ GM*M&/ R*, require large values of D, that in [T2] 
were associated with a postulated turbulent behavior of 
the disk fluid. Since by definition L = uir 2 , we can infer 
E) = 2Y,v t H/3, vt being a turbulence velocity, expressible 
as v t = av s , where a is a free parameter. The fundamen- 
tal question is whether the axisymmetric disk is linearly 
stable with respect to small perturbations that preserve 
its symmetry. A solution to this problem comes from the 
presence of a non- vanishing magnetic field, which makes 
the non-linear interaction of very small disturbances of 
the equilibrium possible. Such a MHD instability, com- 
monly known as Magneto- Rotational Instability (MM), 
is triggered by the radial gradient of the disk angular 
velocity and has been fixed by the Velikov analysis of 
1959 [HI E] (for a review on the topic, see [IT]). In [20] it 
is argued that the MRI is strongly suppressed when the 
disk is thin enough and therefore its efficiency in generat- 
ing turbulence is ruled out in favor of a thermorotational 
instability, in which the vertical gradient of the temper- 
ature plays a crucial role. 

Indeed, the powerful scenario emerging from the MRI 
relies on the presence of an even small, embedded mag- 
netic field. However, the existence of an unstable mode 
of wavenumber k is subjected to the condition 

,22 da: 2 ^ n . . /GM* 

where va denotes the background Alfven velocity. This 
simple condition can be easily recovered when the space 
dependence of the perturbation is parallel to the direc- 
tion of magnetic field, here assumed without a signifi- 
cant loss of generality, along the z-axis. Condition ^ 
can be recast in a more significant shape, as far as we 
realize that the stability of the disk profile as a whole 
requires the extension of such inequality to scale of the 
disk depth. A global disk stability is, in fact, ensured 
only if the Alfven term dominates even for the smallest 
available wavenumbers k ~ ttu}k/V2v s . Hence, observing 
that dw^/dlnr = — 3w%_, we finally get 

vl < \v 2 s . (3) 

As far as the magnetic field being sufficiently small and 
the disk thick enough, the emergence of an unstable mode 
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is guaranteed. However, if the magnetic field of the cen- 
tral object is important, so that the Alfven velocity is not 
too small and at the same time, the temperature and an- 
gular rotation of the disk are able to produce a sufficiently 
thin profile (for which the sound speed is constrained by 
a small upper bound v B /(ojjgr) ~ H/r -C 1), behaviors 
violating condition ^ must be taken into account. In 
other words, we cannot exclude the existence of a class 
of thin disks for which the plasma admits a parameter j3 
less than a few units. Indeed, in the isothermal case we 
have 



2-§- < — 
vl ~ 3 



(4) 



For such relatively cold and strongly magnetized 
plasma disks, the MRI is suppressed and we need a new 
type of instability to account for the turbulent scenario 
required to deal with significant dissipative effects. In- 
teresting features in favor of new instability perspectives 
are discussed in 

The presence of an intense magnetic field (as required 
by a pulsar-like central object) however, poses a new puz- 
zle regarding the value of the non-zero resistivity coeffi- 
cient rj of the ideal plasma. In fact, the equation of the 
electron force balance reads as 



The puzzle consists of the fact that the strong dipolc 
magnetic field is associated with a zero density current 
(being a vacuum solution), and thus the toroidal cur- 
rent J<£ must be rather small since its existence is due 
to the plasma backreaction only. The key quantity we 
have to focus on is therefore the ratio Mjfio/i), which is 
requested to be very small as well. For an X-ray star, for 
which the parameters and fio, as determined from the 
observations, are significantly high, the smallness of the 
toroidal current can be reached only for a sufficiently 
high value of the resistivity coefficient. But this is clearly 
not the case, as far as we microscopically estimate this 
coefficient, accordingly to the expression (in Gauss CGS 
units, i.e., in seconds) 



(10) 



where m c is the electron mass, v pc is the proton-electron 
collision frequency, n D is the electron number density, 
and e is the electron charge. The proton-electron col- 
lision frequency v pc can be well approximated through 
the electron-electron collision frequency f cc , given bjF] 



E 



A B = r/J , 



(5) 



with E and B denoting the electric and magnetic field, 
respectively, while v is the plasma velocity and J is the 
current density. Since the axial symmetry requires E^ = 
0, the azimuthal component stands as 



2.91 x 10" 



( Uc ) 
\ 1 cm -3 / 



T 
TeV 



log(A)Hz, (II) 



v?B r — v r B 7 



(6) 



In the one-dimensional model, obtained for the thin 
disk by averaging along the z-axis, and taking into ac- 
count both the relation between v r and the constant ac- 
cretion rate , as well as the dipole-like morphology of 
the magnetic field, Eq.Q is given by 



where the Coulomb logarithm log(A) can be estimated 
via the formula 



A = I 2tt7V, 



D 



(12) 



M d /x 
27rS(r)r 4 



(7) 



Ho being a parameter modulating the intensity of the 
dipole field. Observing that the superficial density is 
given by [5] 



N-£, being the Debye number of the plasma. In FigJT] 
this resistivity coefficient is plotted against the number 
density and the temperature of the plasma, showing the 
smallness of the corresponding assumed values. 



E(r) 



2ttv s 2 
GAh 



1/2 



Po(r) 



,3/2 



(8) 



the azimuthal current density takes the radial behav- 
ior 



Jrb 



Md^o / GM» 



1 



rjcv s y (27r) 3 po^r 11 / 2 



(9) 



1 see the NRL Plasma Formulary (Naval Research Laboratory, 
Washington, D.C., 2009). 
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Figure 1. Resistivity coefficient rj as a function of the elec- 
tron number density n e (cm -3 ) and of the temperature of the 
plasma T (eV). 

The viscoresistive puzzle is therefore the combination 
of two related facts: a magnetic field is required by the 
MRI, which allows the onset of turbulence, to generate 
the viscosity term (without such dissipative effect no an- 
gular momentum transport is permitted); this magnetic 
field in turn requires an anomalous resistivity to get a 
satisfactory balance for the electron force along the az- 
imuthal component. We are led, thus, to postulate both 
a high viscosity and resistivity coefficients, i.e., a mag- 
netic Prandtl number of order unity [7J. However, no 
firm explanation exists for the possibility that the turbu- 
lent behavior, resulting from the MRI, is able to generate 
significant dissipative effects in a thin disk, and, as dis- 
cussed above, some plasma configurations can exist for 
which the MRI does not work at all. 



B. The crystalline structure: from the local toward 
a global model 

This apparent inconsistency of the standard model for 
stellar accretion (at least for sufficiently magnetized as- 
trophysical sources) and the possible failure of the MRI 
for a sufficiently thin disk, led Coppi in [TJ [5] to reformu- 
late the thin disk equilibrium in a pure two-dimensional 
ideal MHD scenario. The point of view proposed cor- 
responds to replacing the diffusive morphology of the 
magnetic field in the plasma, resulting from the effec- 
tive dissipative processes, by rigid structures in the mag- 
netic flux surfaces and in the mass density profile, the 
so-called crystalline structure, with the associated ring 
sequence decomposition of the disk. In this scenario, the 
angular momentum transport is not ensured by the tur- 
bulent viscosity effect, but is argued to come from the 
porosity of the plasma near the A-points of the mag- 



netic configuration, where the z-component of the B field 
vanishes. Indeed, near such peculiar configurations, the 
radial velocity can increase without any violation of the 
electron force balance. In such a framework, the plasma 
is pushed to infall by the action of intermittent thermoro- 
tational unstable modes, similar to the ballooning modes 
observed in laboratory experiments [231 24 ; for a dis- 
cussion of this inferred alternative scenario for the disk 
accretion, see [It)] . 

We also stress how the possibility to reconcile the crys- 
talline structure with the viscoresistive MHD scenario is 
discussed in [35], showing how the radial oscillation of 
the magnetic flux surfaces can be recognized, for a local 
model, even in the presence of dissipative effects. Fur- 
thermore, [26 studied the ideal MHD configuration of a 
thin disk near a fixed radius from the central object, but 
differently from the original works [HE], the contribution 
of small poloidal currents and matter fluxes are retained 
in the problem. The idea proposed concerns the possi- 
bility to balance the electron force equation by requiring 
the peculiar condition v z B r — v r B z ~ (to be intended 
as a net average prescription). In the framework of a 
plasma decomposed in a ring sequence, the possibility to 
deal with strong values of the vertical velocity, in corre- 
spondence to suitable restrictions, is shown. The estab- 
lishment of such configurations seems to be of relevance 
in the explanation of stellar wind profiles, as well as in 
the determination of seeds for the jet formation. 

All these results on the crystalline structure, however, 
are derived within the scheme of a local model, in which 
the plasma configuration is fixed around a fiducial radius 
from the central object. The main goal of our analysis 
is the extension of the original local model to the global 
configuration of the disk. 

Indeed, the local analysis is developed around a fixed 
value of the radial coordinate in the disk, and therefore it 
is unable to clarify how the differentially rotating layers 
determine the radial behavior of the fundamental physi- 
cal quantities, such as the pressure, the mass density and 
the magnetic flux surfaces. One of the main questions left 
open is whether the crystalline structure and the ring se- 
quence arising in the disk, are compatible with a regular 
radial envelop, i.e., if these peculiar features can be rec- 
onciled with a physically meaningful behavior of the fun- 
damental configuration variables. For instance, an im- 
portant point to be addressed, and successfully handled 
in this paper, is the possibility to have a radial oscilla- 
tion of the mass density (characterizing the ring forma- 
tion) but in the framework of a global decay of the radial 
matter distribution as the outer regions of the disk are 
approached. 

Having this scenario in mind, we formulate the global 
disk plasma equilibrium starting from those natural as- 
sumptions which are the expected generalization of what 
was proposed in [HE]. As first step, we need to fix the 
explicit dependence that the plasma angular frequency 
outlines, according to the corotation theorem [27], with 
respect to the background magnetic surface; the latter, 
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in particular, are taken in the suitable form of a dipolar- 
like configuration. In the local model, this information 
is somehow hidden in a constant term and ultimately 
restated as the derivative of the angular frequency with 
respect to the radial variable evaluated at the fixed point. 
Here we fix the dependence w = uj(tp) by determining it 
on the equatorial plane of the disk and then extending 
this relation to the whole plasma profile. 

Another basic assumption at the ground of the global 
analysis presented below, is the request of dealing with 
a very small-scale backreaction of the plasma with re- 
spect to the magnetic field of the central object. In other 
words, we assume that the physical quantities (i.e., the 
internal currents, the mass density, pressure perturba- 
tions and the correction to the magnetic surfaces) have a 
very rapid radial variation, in agreement with the MHD. 
This request is equivalent to fixing a hierarchy in the 
radial gradient profile, and has the same impact on the 
plasma configuration as the smoothness of the plasma 
around the fixed radius, which is implicitly postulated in 

These two requirements at the ground of the global 
model, i.e., the dipole like nature of the central object 
magnetic field and the very small scale of the plasma 
backreaction, appear as very natural choices for the phys- 
ical characterization of a stellar disk, and their success- 
ful implementation to obtain the crystalline array of the 
plasma along the radial direction, must be accompanied 
only by the restriction to deal with a thin disk; such an 
assumption (already present in the local analysis) is al- 
lowed by the many observations that such a class of disks 
is indeed present around compact and rapidly rotating 
astrophysical objects. 

The real new contribution of what is discussed here, 
is the determination of the radial profile of the disk and 
the estimation of the microscopic nature (for an astro- 
physical setting) of the wavelength of the perturbations 
responsible for the ring decomposition of the disk. Such 
an analysis has no direct relation with [25], where the 
poloidal currents and matter fluxes are taken into ac- 
count; however, the natural development of the present 
global model would be to verify the existence of such 
peaks in the vertical velocity, once the radial envelope is 
fully accounted for. This scenario would allow the ex- 
act determination of the disk regions where the stellar 
wind or the jets could take place. Nonetheless, involv- 
ing poloidal currents and matter fluxes in a global model 
would lead to non-trivial questions concerning the valid- 
ity of the corotation theorem, well beyond the scope of 
the present approach. In fact, in the addressed toroidal 
picture, the validity of such a theorem stands as well and 
we are focusing our attention on the radial dependence of 
the crystalline structure, i.e., on how the magnetic flux 
surfaces and the mass density oscillate across the global 
radial profile of the disk. 



III. FUNDAMENTAL EQUATIONS 

We now fix the fundamental configuration equations 
for the axial symmetry of the disk, by considering the 
central astrophysical object as a star of mass M*, en- 
dowed with an intense magnetic field B. The magnetic 
field can be taken in the form 



B 



1 



d z i>e r 



1 



d r ipe z 



(13) 



ip = ip(r, z 2 ) being the magnetic flux function, while the 
associated current density J remains fixed as 



J = 



c 

47T 



8 r 



1 



(14) 



The magnetic field of the central object is well described 
in the disk by a dipole-like configuration, corresponding 
to the flux function 



ifo{r, z 2 ) 



( r 2 + z 2)3/2 



Ho = const . (15) 



Here, the constant value /io fixes the dipole field ampli- 
tude and we stress that the current density associated 
with this vacuum configuration is characterized by the 
relation J(ip = ?/>rj) = 0. 

The Newton potential x describing the gravitational 
field generated by the central object stands as 



X(r, z 2 



GM, 



\Jr 2 + z 2 



(16) 



G being the Newton constant. We can also define the 
Keplerian angular velocity u>k as 



w|(r, z 2 ) 



GM, 



(r 2 + z 2 ) 



2^3/2 



(17) 



We observe that on the equatorial plane z = 0, the fol- 
lowing relations hold 



w|(r,0) = w| = 



GAL 



2 _ GM, 



(18) 



These relations are of interest because of the thin na- 
ture of the disk. Furthermore the corotation theorem 
[27] states that the plasma angular velocity must be a 
function of the magnetic surface only, i.e., u = u>(t/j), 
and it is natural to postulate that the relation 



GM* 
Mo 



(19) 



holding on the equatorial plane, is valid everywhere in 



the disk. We note that, in agreement with Eq.(17l, the 



disk, embedded in the dipole magnetic field of the central 
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object, can not have a Keplerian behavior far from the 
plane z = 0. 

The radial equilibrium of the disk configuration corre- 
sponds to the following force balance 



GM,p 



rip 3 
Mo 



ipr> 
Por 



= -d r p 



d r ip 
Anr 



Vty] , (20) 



where p — p(r, z 2 ) and p — p(r,z 2 ) denote the mass 
density and the thermostatic pressure, respectively. 
The vertical configuration equation reads as 



d z p- p 



GM^z 
p r 2 



d z ip 
4irr 



vm = o . 



where 



= d r 



-d r ip) + -d 2 ^. 

r I r 



(21) 



(22) 



IV. PERTURBATION SCHEME 



We now split the flux surface function ip as follows: 



(23) 



Here £(r, z 2 ) describes the electromagnetic backreaction 
of the confined plasma. Via the approximation tp 3 ~ 
■0q + 3-0^7 an d recalling that the toroidal currents, asso- 
ciated with ip-Q identically vanish, the radial and vertical 
configuration equations rewrite as 



GM, 



ipu 



--2(^ + 3^C) 
Mo L Mo 

d rP - ^:(d r tp B + d r C)V[C}, 



(24a) 



Airr 



, , GM,z 

OzP + P ^~1pD 



Anr 



{d z ifo + d z Q V[C] = , (24b) 



respectively. In these equations, we retained the gradi- 
ents of C, because, in agreement with the drift order- 
ing, they can be relevant and even dominant, despite the 
smallness of £. 

Taking into account the thin nature of the disk, i.e., 
that its half-depth H (r) satisfies the relation H(r) <C r, 
we can expand the vertical dependence of the above 
equations in agreement with the approximation (r 2 + 
z 2 ) a ~ r 2a (l + az 2 /r 2 ). Hence, the configurational sys- 
tem Eq. (24 1 is recast as 



ZpGhh 



(r 3 - 3rz 2 )( 



Mo 



-d r p 



1 

Aur 



d z p + p 



GNhz 



Mo 



1 

4irr 



_ 9^ 
1 ~ 2T 2 " 

3^o^ 



(25a) 



V[(] 



d z {)V[t]=0. (25b) 



We now split the mass density p and the pressure p into 
the background (overbar) and perturbation (careted) 
components, as follows 



P = P + P, 



p = p + p, 



(26) 



requiring that the quantities p and p are linked via the 
isothermal relation p = v 2 p. We also determine the form 
of p(r, z 2 ) by imposing the validity of the gravothermal 
vertical equilibrium 



d z p + P^kqZ = 



p (r)exp - 



2v 2 



(27) 



where we remind that po(r) denotes the mass density on 
the equatorial plane. Observing that 



d r p 



(ZGM*z 2 
\ 2r 4 



v^dpo 
Po dr 



(28) 



in order to restate the radial and the vertical equilibria 
in a simpler form, we are led to require that the following 
conditions hold 



Mo 



C 3Z 2 

0do >> 2r 2 ' 



« 1, 



C rv 2 

— ^— > 7 — 

V>do 3GA7* 



(29) 



("0DO = ip~D(r,z = 0)). Here we make the ansatz po oc 
r 7 , 7 = const (see below) and obtain the following system 



of equations from Eqs.(25) 



(P + P) 



3GAL 



= r p- 



1 



d z p + p 



GM, 



Mo 
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Airr 
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1-**) 

2r 2 



3poz 



1-3- 



(30a) 



d z ()V{(}=0. (30b) 



A. Toward the crystalline structure 

In order to find the region of applicability for the crys- 
talline structure of the disk, outlined in the local model 
of Coppi Q], we impose some specific restrictions. De- 
noting by k the wavenumber of the radial dependence 
characterizing cp = C/V^i we make the request to live 
in the disk zone where > 1. Thus, under such a 
restriction and neglecting the quantity d z ipu in the ver- 
tical force balance, the disk configuration is determined 
via the following system 



3GM* Po (r) 



PoVr 



exp 



- Ku - 



2v 2 



d r p + 



4nr 



Mo 

r 5/2 



D 



2r 2 ) 



1 



{d 2 r <P + d 2 z <p) 



(31a) 
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d z p- 



GAL 



+ ^d z( f>(d^ + d^)=0, (31b) 



where we defined D = D(r,z 2 ) = p/po- If we indi- 
cate by h the characteristic scale for the z-dependence 
of C, the possibility to neglect the term d z ipD in the ver- 
tical equilibrium, relies on the validity of the requirement 
C/V-D > izh/r 2 . 



V. THE LINEAR CASE 

We now study the linear case, corresponding to a suf- 
ficiently low electromagnetic backreaction in the plasma, 
so that we can neglect non-linear terms in ip when fixing 
the equilibrium. For instance, d r ipu 3> d r (, is equivalent 
to the condition 



■0D ^ kr 



(32) 



Furthermore, we observe that the background mass den- 
sity can be expanded as 



Po(r) 1 



GALz 2 
2r 3 v 2 



Po(r) 



Rsc 2 
4rv 2 



v2 1 



(33) 

where i?s = 2GM*/c 2 is the Schrwarzshild radius of the 
central object. In what follows, we assume to be in that 
region of the disk where the condition 



Rs v. 
— > - 
4r c 



(34) 



holds and we can approximate, in this linear regime, the 
configurational system (31) as 



SGAL 



Po- 
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L s z 2 



d z p + p 



Po 
4 7rr 3 



GAL 



(d 2 ct> + d 2 (j)) =0, (35a) 



= 0, 



where 



L s ee GM*/2v, 



2 

s ' 



(35b) 



(36) 



(see Fig|2J) . In these equations we neglected D, because 
p ^> p in the linear regime, and also the radial pressure 
gradient, in comparison to the Lorentz force. We observe 
that condition (34) reads as r <C L s and that it is im- 



plied by condition ( 29 ) . Equation ( 35b ) tells us that the 



perturbations follow t 
corresponding radial c 
feature: if we take po 



;he background behavior, while the 
corresponding radial equilibrium ( 35a ) has an intriguing 

m/r 3 , 



with 



= const, and 




T(K) 



Figure 2. The parameter L a as a function of the temperature 
T of the isothermal disk for a M, = 2M Q 



define k 2 = Yl-KGM^mj 'pq, then we arrive at the radial 
equilibrium equation in the form 



(d 2 4> + d 2 cb) 



k 2 (j). 



(37) 



In the considered limits z/r <C 1, kr 3> 1 and for r 3> 
(-L s //c 2 ) 1//3 (and remembering r <C L s ), the solution of 
the equation above takes the form 



r ,z 2 ) = A sin(fcr) exp 



k^/T s z 2 
2r 3 / 2 



(38) 



A being a constant amplitude. This oscillating form of 
</>, i.e., of C = \/rcj) restores, under the set of conditions 
we fixed, the crystalline structure emerging in the local 
model of Coppi. The local oscillating behavior of £ is 
sketched in Fig(3j while the global behavior is described 
in Fig|4j 




Figure 3. Local oscillating behavior of the perturbed flux 
function £ around a radius r = 10 3 km. The chosen parame- 
ters are: A = 1, T = 10 7 K, L B ~ 8 x 10 5 Km, B = 10 12 Gauss 
and m = 0.02M Q 

All the conditions imposed on the global model to 
recover the linearized radial oscillating behavior of the 
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the wavelength A = 2n/k and the typical scale L s , i.e. 



2xl0 6 5xl0 6 lxlO 7 2xl0 7 5xl0 7 



R(m) 



Figure 4. Global profile of the perturbed flux function £ for 
fixed values of z = 0, 10 2 and 10 3 m (clearly, the oscillating 
factor sin(fcr) has been omitted). The chosen parameters are 
the same as in Fig[3] The weak radial dependence of the 
amplitude is evident. 



magnetic surface function, can be summarized in the fol- 
lowing three physical restrictions on the considered disk 
model: 

1. The disk is assumed to be thin, i.e., z 2 /r 2 <C 1. 

2. The wavelength of the perturbations, due to the 
electromagnetic backreaction, is much smaller than 
the disk size, i.e., kr ^> 1. 

3. The disk configuration must be restricted to the 
radial region where the condition r <C L s holds, 
i.e., 



r < 



G¥» R s c 2 



2v 2 



4v 2 



m x cr 
8K B T 2 



Rs 



(39) 



where, the sound velocity v s has ben expressed via 
the disk temperature T. 

In this scheme, the existence of the crystalline struc- 
ture requires that the disk temperature not exceed a 
given value. For instance, for a neutron star of about 
two Solar masses and a disk size of 10 6 km, in order to 
apply our model over the whole disk, we would get the 
condition K B T <C IkeV. Such an estimation leads us 
to infer that the present result, as that obtained in [T], 
remains valid for a class of rather cold disks. 



VI. PHENOMENOLOGICAL IMPLICATIONS 

Let us analyze the physical hints that we can get from 
the analysis of the global model for the crystalline struc- 
ture, as traced above. 

In order to compare the radial scaling of the pertur- 
bations to the size of the allowed region for the crys- 
talline structure, we have to consider the ratio between 



A = / H 

L s V 3GM*m GAL 




-3-^/ off. 



where we defined the characteristic length 



(41) 



Since v s < c and t -C Rs (for a typical neutron star 
we have value in the range ~ 10m; see below for more 
details), the above ratio is many orders of magnitude 
smaller than unity. The oscillation scale of the pertur- 
bations, therefore, lives in the microscales of the system, 
very much below the characteristic length L s . The most 
relevant phenomenological implication of this issue, is in 
the impact that the crystalline structure can have on the 
global configuration of the disk. The very small char- 
acteristic length of the magnetic flux surface oscillations 
prevents the ring sequence emerging in the non-linear 
regime (see Sec. VII) from being on a macroscopic scale, 



and so able to account for certain discontinuous struc- 
tures resulting from the observations of real astrophysical 
compact objects. Such kind of macroscopic structures 
can take place in this framework in limiting situations 
only, i.e., when the parameter m is particularly small 
and fiQ is sufficiently large. However, this consideration 
does not affect the relevance of the existence of such mi- 
crostructures toward the global equilibrium profile. In 
fact, this microscopic nature of the crystalline structure 
and of the ring sequence, plays a crucial role in fixing the 
fundamental instabilities characterizing the disk profile. 
In particular, the small-scale radial oscillations can trig- 
ger instabilities having very different morphologies with 
respect to the MRI, as already discussed for the local 
configuration in |20j . Thus, more than through direct 
observations in the optical or X-ray bands, we expect 
that the existence of such structures inside the disk can 
be revealed by specific features of the triggered turbu- 
lent regimes that can be unveiled in the resulting emis- 
sion processes. The role that the microscopic crystalline 
structure can play in the onset and in the establish- 
ment of a turbulent scenario, can deeply influence the 
transport processes inside the disk plasma, especially in 
view of energy-momentum transfer from the micro- to 
the meso-scales of the system. 

Another relevant phenomenological question concerns 
the restricted region where the crystalline profile and the 
fragmentation of the disk in a ring scries can take place. 
Indeed, the condition r <C L s seems to confine the ra- 
dial domain allowed for the magnetic flux surface oscil- 
lations in an internal portion of the thin disk. From a 
phenomenological point of view, this inequality is sensi- 
tive to the accretion system parameters, i.e., the mass 
of the central object and the temperature of the disk, so 
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identifying which fraction of the plasma extension [^] is 
involved in the crystalline profile. In particular, for con- 
figurations extended enough, having very large values of 
i? cxt and whose plasma is significantly hot (so that the 
sound velocity is a few orders of magnitude smaller than 
the speed of light), the restriction above can limit the 
small scale radial oscillations to a very tiny portion of 
the accreting plasma. However, it is possible to charac- 
terize the class of disks to which such restrictions can be 
applied in a more simple and meaningful way, which is 
also consistent with the approximation scheme adopted 
above. In fact, the restriction r <C L s can be easily re- 
stated as follows 



r < L s = 



GM, 
2vf 



2v* 



r 

HHr) 



(42) 



where we recall that H(r) is the effective half-depth of 
the isothermal disk (see Eq.(27l). Thus, the restriction 
r <C L s is equivalent to the requirement that the disk 
has a real thin profile, i.e., H(r) <C r. Realizing such 
an equivalence allows us to join together two of the fun- 
damental constraints at the ground of our derivation of 
the global model, i.e., r -C L s and z/r -C 1. In fact 
the former, being equivalent to H(r) <C r automatically 
ensures the validity of the latter. It is worth noting, how- 
ever, that for a thick disk, while the constraint r -C L s 
is necessarily violated (with non-trivial implications for 
the existence of a global crystalline profile), the condi- 
tion z <C r can still be satisfied in the proximity of the 
equatorial plane of the axisymmetric configuration. 

We conclude this section by stressing how the two fun- 
damental requirements we relied on when deriving the 
radial oscillation of the disk morphology, overlap the fun- 
damental hypotheses of dealing with a thin disk in which 
the plasma backreaction is a small scale phenomenon, 
namely, H(r) <C r and kr ^> 1, respectively. This issue 
makes our analysis fully consistent with the ideas intro- 
duced in [U |2] and substantiates the guess that the crys- 
talline morphology and the ring sequence scenario are 
very general features of an accreting plasma well con- 
fined close to the equatorial plane. However, our study 
also has the merit to outline the microscopic nature of 
the backreaction scale (more than a simple small scale 
behavior postulated by Coppi), so opening a precise di- 
rection in the understanding of the role that such a radial 
periodicity of the plasma can play in the establishment of 
significant processes of transport of matter and angular 
momentum across the disk. 

With respect to the estimation of the characteristic 
length £, it is worth noting the following relation (valid 
in the linear case) between the total mass of the disk Md 
and the parameter m entering the equatorial distribution 



2 Such fraction is roughly proportional to the ratio L s /_R ex t, where 
i?ext denotes the external radius of the disk configuration. 



p (r) = m/r 3 



Md = 2irm 



H(r) 



H(r) 



1 



e-* 2 ' H2 ^dzdr 



1/2 
int 



(43) 



1.8 x 10 3 m 



kcV 



Here erf (a;) = 2tt~ 1 / 2 J® e _t dt is the error function, 
M* = 2M© is the mass of the central object, and the 
last equality stands for an inner and external disk radii 
of i?i n t = 10 3 km and i? C xt = 10 6 km respectively. As a 
result we find that the parameter m explicitly depends 
on the temperature T and on the mass of the disk Md 
via the relation 



m 



5.5 x 10~ b M d 



K B T 
kcV 



(44) 



Assuming that Md takes values in the range 
[1O- 5 ,1O~ 2 ]M , K B T may vary within [10" 1 , 10 2 ]keV, 
the stellar radius to be equal to 10km, the parameter m 
assumes values in the interval [5.5 x 10 -11 , 1.7x 1O _5 ]M . 
The parameter £ given by Eq.(41|, then, is related to 



the parameters of the model as follows 



2.3 



/ b y ( K B T 
\10 9 GJ \keV 



M d 

M 



(45) 



For different values of the mass of the accretion disk Md 
and for a B-field of 10 12 Gauss, the behavior of £ and the 
ratio A/L s ( |40| as functions of the temperature are given 
in Fig[5]and Fig|6] respectively. From the behavior of I, 
it is evident that the condition I i?g is always verified. 




100 



T (KeV) 



Figure 5. Behavior of the characteristic length I as a function 
of the temperature T of the disk. The different values above 
each single plot denotes the mass of the accretion disk Ma in 
Solar mass units for which Eq.(|45[) is evaluated. It is worth 
noting how such a characteristic length is always much smaller 
than the typical dimensions of a disk. 
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50 100 



Figure 6. Plot of the X/L s ratio as a function of the 
temperature T of the disk. The different curves corre- 
spond to different values of Ma /Mb (fr om top to bottom: 
10" 5 , 1CT 4 , 10" 3 , 10" 2 ) for which Eq.Jiol is evaluated. 



VII. EXTREME NON-LINEAR REGIME 

We now analyze the opposite case with respect to 
the linear regime, when the electromagnetic backreac- 
tion dominates the background dipole contribution (i.e., 
d r Q 3> dr^i)) and the perturbed mass density p is larger 
than p fixed by the gravostatic equilibrium (33); more 



precisely, we can apply the following hierarchy relations 



p~»p. 



(46) 



Therefore, we can also require that the gravitational term 
of the system (31 ) be negligible with respect to the pres- 
sure gradient, i.e., (here we are assuming p = v 2 p) 



r 2hz 

— > — 2" 



(47) 



where L s = GM*/2v 2 . We now split the pressure term 
in analogy to the surface function (, i.e., p(r, z 2 ) = 
q(r, z 2 )/r, so that 



d r p = - (d r q - -q) ~ -d r q , 
r V r I r 



(48) 



valid for kr 1. Hence, fixing the mass density on the 
equatorial plane in the form po = m! /y/r, the configura- 
tion system ([311 reads as 



3GM,m' ^ 1 
D<p = d r q + —a, 

Po 47T 

d z q+^d z <f>A<f> = 0. 

47T 



(49a) 
(49b) 



In order to analyze the system above, let us introduce 
the following dimensionless quantities 



x = kr 



z 
h 



i 



s = 



(kh) 2 
3GikL™ x 2/9 



(50) 



T (KeV) w here we set ml 



/X fc 3 / 2 ' 

m/L 5 / 2 , with L — const. By such 



definitions, system ( 49 1 rewrites in the form 



M + d x s + ~d x $ (a 2 $ + tfl®) = o , 



9uS + — 9 U $ (a 2 $ + ea 2 $) =o. 



(51a) 



(51b) 



We now solve these two equations in close analogy to 
what is done in [2] in this same extreme non-linear 
regime. We consider a function $ in the form $>(x, u 2 ) = 
N(x)F(u 2 ), where N is an odd function of the radial 
variable. In the limit £ — > 0, which we are focusing on, 
from Eq.( 51b[ ), we easily get 



S(x, u 2 ) 



1 2* r d 2 N 

F 2 N — - 

8tt dx 2 



(52) 



Substituting this expression into Eq.(51a) and setting 
D(x,u) = K(x)F(u 2 ), we eventually get 



8irK{x) 



d 3 N 1 dN d 2 N 



dx 3 N dx dx 2 



(53) 



where we must necessarily require K > 0, to ensure the 
positivity of the mass density p = DpQ. 

If we assume the following simple form for N(x) 



then 



N(x) = A [sin(x) + B sm{2x)] , (54) 

, . 6AB sin 2 (x) , , 

8ttK(x) = — , . v ' . (55) 

v ' 2Bcos(x) + l v ; 

Finally, we can calculate the quantity v 2 via the relation 



. 2 _ p _ p 2 k 1Y l 2 ti>l 2 



AF (u 2 ) [2Bcos{x) + l} 2 [8Bcos{x) + 1] 

6By/x~ 



(56) 



Since we have to require K{x) and v 2 to be positive defi- 
nite, then parameters A and B must satisfy the following 
conditions 



A > 0, < B < 



1 



(57) 
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For the particular choice A = 1 and B = 1/8, and 
F(u 2 ) = 1, the (effective) perturbed sound velocity i) 2 
( 56 ) and the perturbed mass density 



P = Dp = 



K(x)fh 



3m sin 2 (a;) 



L 5 /2 A 



V 8ttL 5 / 2 y/xfk [cos(ar) + 4] 



(58) 



are depicted in FigjTJ The radial dependence of i) 2 can be 
interpreted as a corresponding behavior for the temper- 
ature perturbation induced by the plasma back-reaction. 
For sufficiently small scales of the crystalline structure, 
this temperature contribution becomes the dominant one 
and, in this scheme, the disk acquires a non-isothermal 
profile. 

We stress how, even in this global model, the emer- 
gence of a ring-sequence decomposition of the disk takes 
place as far as the induced magnetic field dominates the 
central object one, i.e., when the perturbations have a 
sufficiently small scale such that d r C, 3> d r ipY). This fact 
confirms that a thin disk admits a plasma configuration 
characterized by microstructures which have the form of 
double and opposite current filaments. The instability 
properties of such radial profile of the disk is expected to 
deeply influence the averaged transport features of the 
global structure, especially in view of energy and angu- 
lar momentum transfer from the micro- to the mesoscales 
of the system. 



VIII. CONCLUDING REMARKS 

We considered a two-dimensional axisymmetric ideal 
MHD model to describe the structure of a plasma disk 
surrounding a compact astrophysical object. We ne- 
glected the self-gravity of the disk and treated the cen- 
tral gravitational field as a Newtonian profile induced 
by the mass of the compact body. The magnetic field 
external to the plasma is described by an exact dipole- 
like structure associated with the intrinsic features of the 
central astrophysical source. We then settled the equi- 
librium configuration of the plasma by accounting for 
the internal electromagnetic back-reaction, described via 
static perturbations of the magnetic flux surfaces and 
of the thermodynamical quantities. This way, we fixed 
the global model associated with the analysis pursued 
in pj [2] , where a crystalline profile and a ring sequence 
for the disk emerged as the result of local configurations 
confined around fixed values of the radial coordinate. 
In this respect, we neglected the poloidal component of 
the matter flux, retaining the disk rotation as the only 
dominant effect. Then the drift ordering was taken into 
account when splitting the configuration variables into 
background and perturbation components. The possi- 
bility to link the present global model to the local con- 
figuration is offered by a suitable radial scaling of the 
amplitudes of the oscillating-like behaviors emerging at 




8nLp 
3 m ^fk 




Figure 7. Adimensional sound velocity fij (top) and perturbed 
mass density p (bottom) for the case A = 1 and B = 1/8, and 
F{u 2 ) = 1. 



a fixed radius; furthermore our scheme relies on the small 
scale nature of the perturbation in comparison to the disk 
size. We were able both in the linear and in the extreme 
non-linear cases, discussed above in detail, to recover the 
corresponding local equilibrium profiles. The amplitude 
scaling of the perturbed magnetic flux surfaces, mass den- 
sity and pressure were properly fixed. These results show 
that the crystal structure of the magnetic field and the 
ring-like profile of the mass density are rather general fea- 
tures of the two-dimensional, axisymmetric, ideal MHD 
model. 

The derivation of the radial scaling requires, however, a 
certain number of restrictions on the fundamental quan- 
tities involved in setting the plasma distribution. A key 
request is that the plasma (suffering this decomposition 
in periodic substructures) must be confined within a cer- 
tain region from the central object, much smaller than 
the characteristic length L s defined in Eq.(36). In the 



case of the most common stellar sources, this request sug- 
gests that the crystal morphology of the disk be prefer- 
ably expected in relative low-temperature plasma disks, 
while the case of a very hot disk profile seems to require 
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a very large value of the central mass. For instance, in 
the case of plasma structures surrounding the very mas- 
sive black holes at the center of Active Galactic Nuclei, 
the central mass is typically of the order 10 8 — 10 9 So- 
lar masses, leading to very high values of the length L s . 
However both in the stellar and in the galactic context, 
the internal nature of the crystal profile, located suffi- 
ciently close to the central object, suggests the necessity 
for General Relativistic effects in describing the gravi- 
tational interaction. For a discussion of the role of the 
present configuration of plasma and its instabilities in the 
Active Galactic Nuclei morphology and for the necessity 
to include relativistic corrections, see [18]. Finally we 



stress how in Sec. |VI| we outlined that such restrictions 
on the radial region where the crystalline structure may 
take place is at all equivalent to requesting the thinness 
of the disk. In this sense, it is immediate to recognize the 
class of accretion disks to which the obtained result can 
apply, i.e., all those configurations for which the sound 
velocity is much smaller than the rotational flow of the 
disk. 

This work was developed within the framework of the 
CGW Collaboration (http:// www.cgwcollaboration.it). 
R.B. acknowledges ICRA (International Center for Rela- 
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work. 
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